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We demonstrate that local modulations of magnetic couplings have a profound effect on the tem- 
perature dependence of the relaxation rate of optical magnons in a wide class of antiferromagnets in 
which gapped excitations coexist with acoustic spin waves. In a two-dimensional collinear antifer- 
romagnet with an easy-plane anisotropy, the disorder-induced relaxation rate of the gapped mode, 
Timp ~ To + A (TlnT) 2 , greatly exceeds the magnon-magnon damping, r m _ m w BT 5 , negligible at 
low temperatures. We measure the lifetime of gapped magnons in a prototype XY antiferromagnet 
BaNi2(P04)2 using a high-resolution neutron-resonance spin-echo technique and find experimen- 
tal data in close accord with the theoretical prediction. Similarly strong effects of disorder in the 
three-dimensional case and in noncollinear antiferromagnets are discussed. 

PACS numbers: 75.10.Jm, 75.40.Gb, 78.70.Nx, 75.50.Ee 



Introduction. — The recent development of the neutron- 
resonance spin-echo technique has led to dramatic im- 
provement of the energy resolution in neutron-scattering 
experiments PHI]. When applied to elementary excita- 
tions in magnetic insulators, this technique allows one to 
measure magnon linewidth with the /ieV accuracy com- 
pared to the meV resolution of a typical triple-axis spec- 
trometer. Damping of quasiparticles depends fundamen- 
tally on the strength of their interactions with each other 
and with impurities, information not accessible directly 
by other measurements. Although theoretical studies of 
magnon damping in antiferromagnets (AFs) go back to 
the 1970s [6], a comprehensive comparison between 
theory and experiment is still missing, mainly due to the 
lack of experimental data. 

Magnon-magnon scattering is traditionally viewed as 
the leading source of temperature-dependent magnon re- 
laxation rates in AFs [5J [B]. Another common relax- 
ation mechanism in solids is the lattice disorder, which 
is responsible for a variety of the low-temperature ef- 
fects, such as residual resistivity of metals [7 and finite 
linewidth of antiferromagnetic resonances |8j. However, 
temperature- dependent effects of disorder are usually ne- 
glected because of the higher powers of T in impurity- 
induced relaxation rates compared to leading scattering 
mechanisms and of the presumed dilute concentration 
and weakness of disorder. The closest analogy is the re- 
sistivity of metals, in which the T = term is due to lat- 
tice imperfections and the temperature-dependent part 
is due to quasiparticle scattering. 

In this work, we demonstrate that scattering on the 
spatial modulations of magnetic couplings should com- 
pletely dominate the low-temperature relaxation rate of 
gapped excitations in a wide class of AFs. Such mod- 
ulations, produced by random lattice distortions, yield 
scattering potential for propagating magnons and, at the 
same time, modify locally their interactions. For an illus- 



tration, we consider an example of the two-dimensional 
(2D) easy-plane AF with one acoustic and one gapped 
excitation branch. In addition to potential scattering, 
responsible for a finite damping r cx of optical 
magnons, see Fig. [ija), there exists an impurity-assisted 
temperature-dependent scattering of gapped magnons 
on thermally-excited acoustic spin waves, see Fig. [ijc), 
which yields r imp (T) cx n,T 2 In 2 T. Despite the pre- 
sumed smallness of impurity concentration n i: at low 
temperatures this mechanism dominates over the con- 
ventional magnon-magnon scattering, Fig. [ljb), which 
carries a much higher power of temperature: r mm cxT 5 . 
We have performed resonant neutron spin-echo measure- 
ments with a few /xeV resolution on a high-quality sam- 
ple of BaNi 2 (P0 4 ) 2 , a prototype 2D planar AF [§]. We 
find that the theory describes very well the experimental 
data for the linewidth of optical magnons. Similar dom- 
inance of the impurity-assisted magnon-magnon scatter- 
ing should persist in the 3D AFs and is even more pro- 
nounced in the noncollinear AFs. We propose further 
experimental tests of this mechanism. 

Theory. — We begin with the spin Hamiltonian of a 
collinear AF with an easy-plane anisotropy induced by 




FIG. 1: (Color online) (a)-(c) Diagrams representing im- 
purity, magnon-magnon, and impurity-assisted scattering of 
the optical magnon (solid lines). Dotted lines are acoustic 
magnons. (d) Schematic energy spectrum of the model ([I]). 
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the single-ion term D > 0: 



f3a /3a terms in the magnon-magnon interaction: 



(1) 



Two examples are the nearest-neighbor AFs on square 
and honeycomb lattices. The latter model, with the non- 
frustrating third-neighbor exchange, is relevant to the 
spin-1 antiferromagnet BaNi 2 (P04)2 {2] discussed below. 

As a consequence of broken XY symmetry, excita- 
tion spectrum in the ordered antiferromagnetic state pos- 
sesses acoustic (a) and gapped {(3) magnon branches: 



2m 



(2) 



see Fig. [ljd) for a sketch. Explicit expressions for c, A, 
and m for BaNi 2 (PC>4)2 are provided in [TU] . 

Defects are present in all crystals. While vacancies and 
substitutions may be eliminated in some materials, inho- 
mogeneous lattice distortions remain an intrinsic source 
of disorder, inducing weak random variations 5 J and SD 
of microscopic parameters in the spin Hamiltonian ([lj 
[llj . Both types of randomness have qualitatively the 
same effect on magnon lifetimes. For example, local mod- 
ification of the single-ion anisotropy SD(Sf) 2 generates 
scattering potential for magnons 
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where c k = a k (/3 k ), Uw = SDS (u k + v k )(u k , +v k ,), and 
u k ,v k are the Bogolyubov transformation parameters. 
For optical magnons at k, k' — > 0, the momentum de- 
pendence is not important, U kk > = 0(SD). For bond 
disorder, all expressions are the same with a substitution 
SD — > 5 J and an additional phase factor, which depends 
on bond orientation and disappears after impurity aver- 
aging. 

For the gapped magnons with k — ¥ 0, scattering am- 
plitude in the second Born approximation, Fig. [TV a), av- 
eraged over spatial distribution of impurities is [1^] 



pimp ~ p n T J 2 m ^ max 

i k ~ioan^, — — 



(4) 



where rij is the impurity concentration, U i = 0(6 J, SD) 
is the averaged impurity potential, and w max is the 
magnon bandwidth [10]. Thus, in 2D, conventional im- 
purity scattering results in a finite zero-temperature re- 
laxation rate of the gapped magnons. 

At low temperatures, the principal scattering channel 
for optical magnons is due to collisions with the ther- 
mally excited acoustic spin waves with cq ~ T <C A. All 
other processes are either forbidden kinematically or ex- 
ponentially suppressed. In this case we can consider only 
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where the first and the second row correspond to the con- 
ventional and to the impurity-assisted magnon-magnon 
scattering, respectively, with Ak = k+ q— q' k'. The 
latter is of the same origin as the conventional impurity 
scattering in ([3| since SD and SJ also modify locally in- 
teractions among magnons |10j . In the one- loop approx- 
imation, ^ and ^ yield the self-energies of Figs, [ljb) 
and (c). Applying standard Matsubara technique, relax- 
ation rates can be expressed as 

qq' 

2 N^,S(Ae), (8) 

qq'k' 



T k = TTH, ^ Kq^q' 



where Ae = £k+£ q -£ q '-£k', Ng,^, = n q (l + n^+n k >) - 
riq>n k ', and n q is the Bose factor. 

There are two important differences between r mm and 
pimp.T in and g First, the total momentum is 
not conserved for impurity scattering. This relaxes kine- 
matic constraints of the 4-magnon scattering processes, 
but requires instead integration over the extra indepen- 
dent momentum k'. Second and most crucial, interaction 
vertices VjJJ^, / and V^ q . k / q / show very different long- 
wavelength behavior as q, q' — > 0. We calculate them 
using the approach similar to [SJ |S] , and find that in the 
long- wavelength limit magnon-magnon interaction |5| is 
V^kq'k'q' ex y/qq' , in accordance with the hydrodynamic 
limit fM] ■ However, for the impurity-assisted scatter- 
ing feL interaction is V k ™ k , q , oc l/^/qq 7 . This can be 
understood as a consequence of an effective long-range 
potential for acoustic magnons produced by the gaped 
magnon while in the vicinity of an impurity. 

The leading T-dependence of r k ™ Q and r^"^ 7, can 
be calculated now using ^ and approximating interac- 
tion vertices with their long-wavelength expressions. The 
main contribution to the integrals in ^ and ^ is de- 
termined by acoustic magnons with q, q' ~ T/c. Then, a 
straightforward power counting yields 



1 k^O ~ B 



T 



(9) 



where B ~ ui max |10j . Thus, the inverse lifetime of an 
optical magnon is proportional to T 5 in 2D. A general- 
ization to higher dimensions gives r mm cx T 2D+1 . The 
T 7 -law for the relaxation rate of optical magnons in 3D 
AFs was previously predicted in [TB] . We note that for a 
given model, the effect of magnon-magnon scattering in 
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^ can be calculated using microscopic parameters, thus 
putting strict bounds on its magnitude. 

The same calculation for r^g T proceeds via the fol- 
lowing integral: 



^imp,T 
k->-0 
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where J = J Q A dq with A ~ 7r/a, Ae = cq — cq' — k' 2 /2m, 

and we used the relation between V^.™^ , , i n (6| and J7kk' 
in pb. The naive power counting in ([To]) already gives 
pimpTT ^ wn jj e a more careful consideration shows 
further enhancement of the scattering as the integrals 
formally diverge [logarithmically] in the q — > region, 
demonstrating an important role of the long-wavelength 
magnons in 2D. This divergence is similar to the one in 
the problem of finite T/v ordering temperature in 2D and 
is regularized similarly by introducing low-energy cutoff. 
The cutoff is either due to a 3D-crossover as in the case 
of some cuprates |15) , or a weak in-plane anisotropy that 
induces small gap ujq in the acoustic branch, the case 
directly relevant to the current work \§\ H3] ■ 

Combining Q and ( 10 1 we obtain impurity-induced 
relaxation rate of gapped magnons 



r imp » r + a 
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where both Tq and A are proportional to rii and to the av- 
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erage strength of disorder U i . As a result, the impurity 
scattering leads to a relaxation rate that carries a sig- 
nificantly lower power of temperature than the magnon- 
magnon scattering mechanism. Therefore, despite possi- 
ble smallness of the combined impurity concentration and 
strength, it should dominate not only the T = lifetime 
of the gapped magnon, but also its temperature depen- 
dence in the entire low-temperature regime. A qualita- 
tive prediction of our consideration is that Tq and A in 



( 10 ) should be of the same order since both terms are 



related to disorder. In addition, for samples of the same 
material of different quality, they must scale with the 
amount of structural disorder in a correlated way. 



In the 3D case, impurity-assisted mechanism ( 10 1 gives 



r™ p ' T oc T 9 / 2 , still dominating the 3D magnon-magnon 
relaxation rate cx T 7 discussed above. 

Experiment. — The experimental part of our work is 
devoted to the neutron spin-echo measurements of the 
magnon lifetime in BaNi 2 (P04)2- This material is a lay- 
ered quasi-2D AF with a honeycomb lattice of spin-1 
Ni 2+ ions and Neel temperature Tn ss 25 K. A compre- 
hensive review of the physical properties of BaNi2(P04)2 
is presented in [9j. Its excitation spectrum has an op- 
tical branch with the gap A « 32 K and an acoustic 
mode, as is sketched in Fig. [ljd) . The fit of the magnon 
dispersion yields the following microscopic parameters: 



Ji = 0.38 meV and J3 = 1.52 meV are exchanges between 
first- and third-neighbor spins, and D = 0.32 meV is the 
single-ion anisotropy. The thermodynamic properties of 
BaNi 2 (P0 4 ) 2 follow the 2D behavior down to T < IK 
and a small gap in the acoustic branch, loq ~ 2K, due to 
weak in-plane anisotropy is consistent with the value of 
the ordering temperature [§]• 

The spin-echo experiments were performed on the 
triple-axis spectrometer IN22 (ILL, Grenoble) by using 
ZETA neutron resonance spin-echo option [17 . The inci- 
dent neutron beam was polarized and the scattered beam 
analyzed from (111) reflection of Cu 2 MnAl Heusler al- 
loy focusing devices. We used a fixed-fc/ configuration, 
with kf = 2.662 A" 1 or k f = 1.97 A" 1 . Different rf- 
flipper configurations were used in order to adapt the 
spin-echo time (energy) £nse (&nse = ^Anse) to the 
magnetic excitation lifetimes, typically in the range of 
5 — 50 ps (130 — 13 ^cV). As for any spin-echo experi- 
ment |181 119] , the measurement of the neutron polariza- 
tion (spin-echo amplitude) after the scattering, -P(iNSE), 
provides us with a direct access to the correlation func- 
tion S (q, ^nse) ■ For a spin- wave excitation described by 
a Lorentzian function in energy of half width T, one can 
show that P(ense) = -Po(snse) exp(-T/£ NS E), in which 
the prefactor Pq depends on the spin-echo resolution. 

For our measurements, we have used a 2 cm sin- 
gle crystal of BaNi 2 (P04) 2 oriented with the a* and c* 
reciprocal axes in the scattering plane. The spin-echo 
data were taken at the antiferromagnetic scattering vec- 
tor Qaf = (1, 0, 0) and the energy transfer AE = 3 meV 
corresponding to the bottom of the dispersion curve of 
the gapped mode [S] . In determining the spin-echo ampli- 
tudes, neutron intensities were corrected for the inelastic 
background, measured at the scattering vector Qaf and 
the energy transfer AE = 5 meV. Results of the tem- 
perature dependence of spin-echo amplitudes for several 
representative Ense's are shown in Fig. [2} Solid lines are 
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FIG. 2: (Color online) Temperature dependence of the polar- 
ization (spin-echo amplitude) of the neutron beam P(T) for 
several representative spin-echo energies. 
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the fits of the spin-echo amplitudes with P = Pqc t / £ nse 
usin g rel axation rate in the functional form given by ([9| 
and (111]), r = r mm + T im P, which we discuss next. Using 



the full set of PIT, £nse) data, experimental results for 
T(T) are extracted from the fits of ln(F) vs £nse at fixed 
temperatures. These results are presented in our Fig. [3] 
together with the theoretical fits. 

Comparison. — The relaxation rate approaches the con- 
stant value of To ~ 25 /ieV at T — > 0, in agreement with 
the expectation Q for the gapped mode in 2D. The 
low-T dependence of the relaxation rate is following the 
power law much slower than T 5 . The quality of the free- 
parameter fit of Ar = r(T)— ro with just the T 5 law is not 
satisfactory for either T(T) or P(T)'s in Figs.[3]and[2] and 
the magnitude of AT also requires an unphysically large 
values of the magnon-magnon scattering parameter B in 
([9]), exceeding theoretical estimates roughly tenfold. On 
the other hand, T 2 In 2 T law gives much more satisfac- 
tory fits in the low- and intermediate-T regime up to 12 
K in both T(T) and P(T), shown as a separate fit by the 
dotted line in Fig. [3} The best fit of T(T), given by solid 
line, is the sum of the magnon-magnon and impurity- 
scattering effects from Q and (11), with the magnon- 
magnon and impurity-assisted parameters B = 15 meV 
and A = 90 /ieV, respectively. The same T(T) is used in 
all three curves of P{T) in Fig. [5J the original data from 
which experimental T(T) is extracted. Magnon band- 
width uj rnax = 64 K and the low-energy cutoff cj = 2 K, 
equal to the gap in the acoustic branch, were used. 

Two remarks are in order concerning the role of the 
magnon-magnon relaxation rate used in Fig. [3} First, fits 
of T(T) in Fig. [3] also include a contribution from scat- 
tering off the thermally excited optical magnons, which 




T(K) 

FIG. 3: (Color online) Temperature dependence of the re- 
laxation rate Y of the optical magnon with k ~ in 
BaNi2(P04)2- Full line is the best theoretical fit including all 
contributions with parameters described in the text. Dashed 
and dotted lines indicate separate contributions of magnon- 
magnon and impurity-assisted magnon-magnon scattering. 



is given by r rr = C ( J) e A / T [TO]- Its contribution is 
roughly equal to that of the T 5 -term {9} at T = 16K 
(= A/2), but diminishes faster at lower T. In the fit 
of r(T) we use the value of C = 260 /xeV, about three 
times the theory estimate: C th » 70 jiteV. Second, the 
theoretical estimate of the magnon-magnon interaction 
parameter in T 5 law ^ is B th ta6 meV, again factor 2.5 
smaller than the one used in the fit (-6 = 15 meV). Alto- 
gether, the magnon-magnon contribution to T(T), shown 
by the dashed line and the corresponding color shading 
in Fig. [3j is likely a generous overestimate of its actual 
role in the relaxation. 

Still, the contribution of the impurity-assisted mech- 
anism in r(T) is very strongly pronounced and is not 
explicable by the conventional scattering mechanisms. 
For example, at 12K the impurity scattering accounts for 
at least 2/3 of the temperature-dependent part of T(T). 



The parameter of the impurity- assisted term in ( 11 ) used 
in the fit is A = 90/ieV, which is of the same order with 
the constant impurity term To, meeting our expectations 
outlined above. This is, again, the strong argument that 
both the constant and the T-dependent terms in the re- 
laxation rate must have the same origin, giving further 
support to the consistency of our explanation of the data. 

The values of A and To cannot be determined theo- 
retically as the impurity concentration and strength are, 
generally, unknown. However, another consistency check 
is possible: the ratio of r to a characteristic energy scale 
of the problem, uj max , should give, according to Q, an 
estimate of the cumulative measure of disorder concentra- 
tion and its strength: m nip (5D/D) 2 ~To/u max ~5- 10~ 3 . 
This translates into a reasonable estimate of the disorder 
and its strength in BaNi 2 (P04) 2 : modulation of mag- 
netic couplings is equivalent to half of a percent of sites 
having SD (6 J) of order D (J). The amount of struc- 
tural distortion in BaNi2(P04)2 [2D] is consistent with 
the magnitude of such variations of magnetic couplings, 
given the strong spin-lattice coupling in this material. 

Other systems. — We propose that similar, and even 
stronger, effects of disorder in the relaxation rate must 
be present in the 2D noncollinear AFs, in which magnon- 
magnon interactions acquire the so-called cubic interac- 
tion terms [21] , absent in the collinear AFs considered 
above. The self-energies associated with such interac- 
tion are the same as in Figs, [ljb) and (c), but with 
two intermediate lines instead of three. With the long- 
wavelength behavior of the impurity interaction to follow 
6V 3 (k, q)ocl/ v /g, as in the considered case, a qualitative 
consideration similar to ( 10 ) leads to: 



k->0 
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where A-i 



''imp 



(5D/D) 2 , an even lower power of T. 
Since the canting of spins can be induced by the exter- 
nal field, we propose an experimental investigation of the 
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effect of such a field on the relaxation rate. For the 3D 
noncollinear AFs we predict T im P' T oc T 5 / 2 . 

Recent neutron spin-echo experiment in a Heisenberg- 
like AF MnF2 Q] have reported significant discrepancies 
between measured relaxation rates and predictions of the 
magnon-magnon scattering theory [5j [6] , precisely in the 
regime of low-T and small-k where the theory is assumed 
to be most reliable. Although the current work concerns 
the dynamics of strongly gapped excitations and our re- 
sults are not directly transferable to the case of MnF2, 
we have, nevertheless, presented a general case in which 
the magnon-magnon scattering mechanism is completely 
overshadowed by impurity scattering, thus suggesting a 
similar consideration in other systems. 

Conclusions. — To conclude, we have presented strong 
evidence of the general situation in which temperature- 
dependence of the relaxation rate of a magnetic excita- 
tion is completely dominated by the effects induced by 
simple structural disorder. Our results are strongly sup- 
ported by the available experimental data. Further the- 
oretical and experimental studies are suggested. 

This work was initiated at the Max-Planck Institute for 
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the Advanced Study Group Program on "Unconventional 
Magnetism in High Fields," which we would like to thank 
for hospitality. The work of A. L. C. was supported by 
the DOE under Grant No. DE-FG02-04ER46174. 
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Spin Hamiltonian 

Here we briefly outline basic steps and main results of 
the spin-wave calculations for the energy spectrum and 
the magnon relaxation rates of the J1—J3 antiferromagnet 
on a honeycomb lattice. The harmonic spin-wave analy- 
sis of the nearest-neighbor Heisenberg honeycomb-lattice 
antiferromagnet can be found, for example, in [T]. 

Geometry of exchange bonds of the considered model is 
schematically shown in Fig.|4j The unit cell of the antifer- 
romagnetic structure coincides with the crystal unit cell 
and contains two oppositely aligned spins S%,i and S2,.; 
in positions (0,0) and p — (a/y/3,0). The elementary 
translation vectors are defined as a! = a( v / 3/2, — 1/2) 
and a2 = a(0,l). The lattice constant in BaNi2(P04)2 
is equal to a = 4.81 A. The reciprocal lattice basis is 
bi = 47r/(\/3a)(l,0) and b 2 = 2 7 r/( v / 3a)(l, The 
volume of the Brilouin zone is Vbz = 8n 2 / V3a 2 . 

The spin Hamiltonian includes Heisenberg exchange 
interactions between first- and third-neighbor spins to- 
gether with the single-ion anisotropy: 

"it = J\ • ($2,1 + S 2ii ;_i + S 2 i ; +3 ) 

i 

+ J 3 £S M -(S 2 ,i+2 + S2,i-2 + $2,1-1+3) (13) 



Here S2,i-i denotes spin in the unit cell R; — ai and 
so on. The microscopic parameters for BaNi 2 (P04)2 




FIG. 4: J1-J3 model in a honeycomb lattice. 



(S = I) were determined from the magnon dispersion 
as Ji = 0.38 meV, J 3 = 1.52 meV, and D = 0.34 meV 
[2J. The second- neighbor exchange was estimated to be 
much smaller J2 = 0.05 meV and is neglected in the fol- 
lowing. 

Applying the Holstein-Primakoff transformation for 
two antiferromagnetic sublattices and performing the 
Fourier transformation 



N 1 / 



1/2 X/ 



Ok , 



b - 1 V" »k(R I+ p) b 
° l - N l/2 Z^ e ° k ' 



(14) 



we obtain the harmonic part of the boson Hamiltonian 
S]T[(3J 13 + C)(a k a k + 6 k fe k ) 

k 

D 



Ho 



(15) 



ak^-k^k + +7 ( a k a -k + &k&-k) + h- c 



where we use the shorthand notations 
J13 = Ji + J3 j 

F k = Ji (e lkl + e lk2 + e lk3 ) + J 3 (e l - + e 
with k„ = k • r„ and 



(16) 



zki 



'k.i 



a n n^ a 
'1 = ^(1.0), r 2 , 3 = - 



V3\ 
2 J ' 



""4,5 



^g(l,±>/3), r 6 = -^(-2,0) 



(17) 



Diagonalization of the quadratic form ([15j with the 
help of the canonical Bogolyubov transformation yields 



n 2 = 



^[ £a (k)4« k + £/ ,(k)^ k 



(18) 



where excitation energies are 



£a(k) = Sy/(3J 13 - l^k|)(3Ji 3 + |F k | + 2D) ,(19) 
epQa) = Sy/(3J 13 + |Fk|)(3Ji 3 - |F k | + 2D) . (20) 

The first magnon branch is gapless, £ Q (0) = 0, and 
reaches the maximum value of 



u&ax = Sy/(2J 3 + 4J 1 )(4J 3 + 2Ji + 2D) « 5.9 meV 

(21) 

at k = [q,q] with q 0.25 in the reciprocal lattice units. 
The second branch describes optical magnons with a fi- 
nite energy gap at k = 

A = 2SV3DJ 13 « 2.8 meV . (22) 



The maximum of the optical branch ui/Z,^ is close to (21 1 
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Long-wavelength limit 

In the long-wavelength limit k — > the energy of the 
acoustic branch has linear dispersion e a (k) » c(ka) with 
the spin- wave velocity 



c = Sy 2( J 3 + | Ji)(3Ji 3 + D) = 4.42 mcV . (23) 
For the optical branch one finds 

{kaf (3#J 13 ) 1/2 



E/»(k)« A + - 



2m 



-S(^3+lJl)(3Jl3- 



, (24) 



with r7i = 0.16 meV" 1 for BaNi 2 (P0 4 ) 2 . 

For small k the Bogolyubov transformation can be 
written explicitly in the following way. First, we trans- 
form from the original Holstcin-Primakoff bosons ai and 
6, to their linear combinations: 



1 



V2 



(dj - hi) , bi 



1 



V2 



(25) 



The Fourier transformed Hamiltonian ( 15 ) takes the fol 
lowing form 

S^[(3Ji3 + £)(4a k + &lA) 

k 



U-2 



(D + |F k |)a k 3_ k +-(D- 



J (26) 
F k |)& k &_ k +h.c. 



2" 1 17 2 

Second, the standard u—v transformation is applied sepa- 
rately for a k and 6 k bosons. In particular, for the acoustic 
branch, a k = UkCtu + Vk<X_ k , we obtain 



1 + d 

2cka 



1 



2cfcc 



1 



cka 
2(1 + d) 

cka 



(27) 



2(1 + d) 

where d = D/(3 J± 3 ) and c = c/(3Ji3). In the case of 
BaNi2(P04)2 the two dimensionless constants are d w 
0.06 and 5 « 0.77. 

Similarly, for optical magnons with k 
K = Uk^k + w k /3l k with 



1 



2^ 



2u 7j 



1 



2sfd 



(28) 



Magnon-magnon interaction 

For a collinear antiferromagnet the interaction between 
spin-waves is described by four-magnon terms in the 
bosonic Hamiltonian. The four-magnon terms of the ex- 
change origin are expressed as 



H 



>x) 



~ J2 F^-*Abtb 2 ai (29) 
1+2=3+4 

— ^ F kl (a\a 3 a 2 bi + b\b 3 b 2 a x + h. 



1+2+3=4 



where a\ stands for a kl etc. The single-ion anisotropy 
contributes 



H 



(an) 



D 

"2iV 



(44 a 2ai + b\b\b 2 b x ) (30) 



1+2=3+4 



D 

4N 



^ {a\a 3 a 2 ai + b\b 3 b 2 b 1 + h. c). 



1+2+3=4 

Performing transformation from a k , fok to a k , /3 k we 
obtain various magnon-magnon terms. The scattering 
of optical (/3) magnons on each other, which will be re- 
ferred to as the roton-roton interaction, can be straight- 
forwardly written as 



V 1 



3 J 13 + D 
N 



1+2=3+4 



(31) 



Derivation of the roton-phonon interaction (scattering of 
the optical magnon on the acoustic one, /3 on a) is more 
involved and we obtain an estimate as 



3(J 3 + |Ji)(l + d) 2 



^2 V k 2fa04a 3 a> 2 /3i 



1+2=3+4 



(32) 



the 



magnon-magnon m- 



(29) and (30) applying the 



The individual terms 
teraction obtained from 
Bogolyubov transformation are proportional to it q M q ' 
and diverge for scattering processes involving acoustic 
magnons, see (27). However, the leading ~ 0{\/^/qq F ) 

0(1) cancel out in their 
sfqq 1 in agreement with 



and the sublcading singularity ~ 
net contribution and V^q. k , , cx 



the hydrodynamic approach [3j 0] . 

Local modulation of magnetic coupling constants due 
to structural disorder, etc., will result in independent 
variations of J- and .D-terms in magnon-magnon inter- 
we obtain action in (29) and (30). Thus, the resultant impurity- 



assisted magnon-magnon interaction will retain the same 
structure as the magnon-magnon interaction, with two 
important differences. First, the momentum in such a 
scattering is not conserved, and, second, the variation 



of J [5 J) is associated only with ( 29 1 and the varia- 
tion 6D will contain only (30) part. Since such vari- 



ations are independent, it suffices to consider one of 
them and treat the associated constant as a free pa- 
rameter. The most important consequence of this con- 
sideration is that, in the impurity scattering, there is 
no cancellation of the individual terms that are pro- 
portional to Uqiiqi cx 1 / y/qq 7 , compared to the case of 



magnon-magnon scattering in ( 32 ) discussed above where 
such a cancellation does take place. Thus, in the long- 
wavelength limit, V^Ei „/ oc l/\/qq', with a coefficient 



T/imp 
V k q;k'q' 

proportional to the impurity concentration and strength 
of the disorder. 
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Relaxation rate of optical magnons 

The lowest-order diagram for the magnon self-energy 
calculated using Matsubara technique is 



£(k, ioj) = 



1 



|Vk(q,q') 



2 ^— ' ico 

q,q' 



(33) 



<^q' ^k' 

x [n q (n q < + l)(nk' + 1) - (n q + l)?v n k'] 
with k' = k + q — q'. Then the damping rate is 

Tk = | ^2 l^(q, q')| 2 [«q( 1 + "q'+ n k') -riq'^k'] 

q,q' 

x <5(e k + e q -e q / - e k ') • (34) 
First, we consider the roton-roton scattering processes. 



The low-temperature asymptote of (34) in this case is 



obtained by taking T <C A and keeping the leading ex- 
ponentially small occupation factor. Then, for an optical 
magnon with k = in two dimensions 



3V 2 

pir . rr 

~ 64^ 



x 6 



/>oo />q p2-7T 

/ qdq q' dq 1 I dip n q 
Jo Jo Jo 



(35) 



VI2 



2m 2m 



2m 



Performing integration in ( 35 1 and using parameters 



for BaNi2(P04)2 discussed above we obtain 



r 1 

J- n 



3(mv; r ) 2 r _ A/T 

64tt 6 



0.035 -r- e 
A 



A/T [meV] . (36) 



Without going into details, there exist another channel of 
scattering that corresponds to a conversion of two rotons 



into two high-energy phonons, /3/3 —¥ aa, which leads to 
the decay rate of the optical mode of the same exponen- 



tial form as in ( 36 1 with a numerical coefficient of the 
same order. 

Finally, the the low-temperature asymptote of the 



roton-phonon scattering in (34) is 



3y 2 r°° pi p 27r 

C = J q 1 2d 1 J q ' 2dq ' J 

x sfcq-cq'- 1 *-^-), (37) 
\ 2m I 



where V lp = V rp \/qq' . Subsequent integration yields 



_3 v 2 T 5 



T 



o - Yo — ~0.18^J [meV]. (38) 
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